Microscopic theory of network glasses 
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Abstract 
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A molecular theory of the glass transition of network forming liquids is developed using a com- 

\ 

bination of self-consistent phonon and liquid state approaches. Both the dynamical transition and 

O ■ 
O ■ 

the entropy crisis characteristic of random first order transitions are mapped out as a function of 
the degree of bonding and the density. Using a scaling relation for a soft-core model to crudely 
translate the densities into temperatures, the theory predicts that the ratio of the dynamical transi- 
tion temperature to the laboratory transition temperature rises as the degree of bonding increases, 
while the Kauzmann temperature falls relative to the laboratory transition. These results indicate 
why highly coordinated liquids should be "strong" while van der Waals liquids without coordination 
are " fragile" . 

PACS numbers: 61.43.Fs,64.70.Pf,65.60.+a 
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While much of the phenomenology of the glass transition is universal, shared by all dense 
liquids, everyday glasses are formed from liquids that possess a network structure at the 
molecular level. The universal features of the liquid-glass transition such as the Kauzmann 



entropy crisis and the Vogel-Fulcher rate Jaw 
underlying random first-order transition 
replica symmetry breaking 



e law 



rave been explained using theories based on an 
like that of mean field models with one step 

3] 



41. Furthermore, the variations of the mean relaxation time 

Q 

and degree of nonexponentiality 3] of relaxation from substance to substance have been 
quantitatively predicted using this approach, provided the configurational heat capacity is 
known as input, thus quantitatively explaining the patterns of behavior verbally character- 
ized using the expression "liquid fragility" p . Detailed microscopic calculations buttressing 
these developments are based on a model of a dense interacting atomic fluid in which the 
particles interact via van der Waals forces alone. Is this starting point appropriate for the 
network glasses? To address this issue, we explore the glass transition in a model of a net- 
work forming liquid. The calculations presented are similar in spirit to those made earlier for 
atomic fluids. Our model assumes cross-linking in a liquid is fluctuating, not quenched as in 
models of rubber [7], taking the mean degree of crosslinking as given from experiment or from 
a separate calculation explicitly accounting for the equilibria between the various different 
bonded species. These calculations suggest a connection between the glassy behavior of 
network systems as viewed by random first order transition theory and rigidity percolation. 
Rigidity percolation features prominently in Phillips's discussion of the role of bonding con- 
straints in glass formation 8|, as well as in the extensive studies of Thorpe^ Q|. A major 
result of the present calculations is that the density of non-bonded interactions needed to 
trigger energy landscape dominated glassy behavior decreases when bonding increases. It 
follows that, consistent with experiment, highly cross-linked network glasses should show 



activated dynamics even at very high temperatures. The theory still predicts a crossover 
to non-activated dynamics for sufficiently high temperature as has recently been inferred in 



computer simulations of silica 
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121. 



The se 
ing point 
functionals 



f-consistent phonon method applied to the amorphous state provides our start- 
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14j . This scheme is closely connected to approaches based on de nsity 
15| . It has recently been justified by an elegant clone or replica approach jrol. 
We actually use a hybrid of both density functional and self-consistent phonon calculations. 
For the hard sphere system, both approaches agree there is a critical density pa, below which 
localized density distributions corresponding to amorphous packings are mechanically unsta- 
ble to even small thermal fluctuations, p^cr 3 from the self-consistent theory is approximately 
1.0 for hard spheres, where a is the hard sphere diameter. This is in agreement with the 



computer simulated dynamical transition 



la] . This instability, from the vitrified side, co- 



incides with the transition found in dynamical mode coupling approaches which start from 



the liquid side |2, 



19 



20]. The thermal vibration allowed at pa in a frozen amorphous 
state gives mean square fluctuations < 5r 2 >= a~ l , where a is an effective spring constant 
and is finite at the transition, a provides an analog of the Edwards- Anderson order param- 
eter for spins and exhibits discontinuous replica symmetry breaking. This limit of stability 
for an amorphous system is comparable to the Lindemann criterion that describes melting 
of crystalline solids. At densities above pa many different mechanically stable amorphous 
packings or free energy minima can exist. At densities beyond pa, however, the configu- 
rational entropy of these minima decreases until it vanishes at a density px- Analogously, 
for soft-sphere models at fixed densities one can find a temperature Ta, signaling the on- 



set of glassy dynamics ("landscape determined" 



21|)and a still lower temperature Tk, the 



Kauzmann temperature where the configurational entropy would vanish. Between these 



densities or temperatures, the dynamics can be described using the concept of entropic 



droplets 



2, 



5[. The activation free energy diverges as T% is approached, consistent with 
the Vogel-Fulcher law AF+ = DT K /(T — T K ). A microscopic calculation of the energy cost 
of forming an entropic droplet quantitatively explains the experimentally observed variation 
of D with the heat capacity drop upon vitrification, AC P 3]. This correlation relies on the 
idea that the Lindemann ratio near Tk should depend only weakly on the detailed force law. 
We will check this supposition for our network glass model. 

The free energy is computed for a given amorphous packing as a function of a. This free 
energy should then be averaged over packings. Each atom can be assigned to a cell and the 
density can be written as p(r) = ^ pi = £\ (^) 3 ^ 2 exp [— a»(r — Rj) 2 ] where the locations 
{Rj}, characterize the packing as a set of fiducial locations around which the atoms vibrate. 
The free energy as a function of {a,}, Fdatj}), is analogous to the inter replica potential 



given as a function of q for spin glasses 



221. 



More generally, the thermal vibrations are described by a tensor a, defined by 



a,. 1 = ((ri-ROfa-RO) 



eff 



(1) 



where the averaging is with respect to the effective potential (dependent on {%}) 



exp [-(3Wf f ) = n^exp (-/?V e// (|r, - Rj|)) = J dv 3 Pj ( Tj ) exp (-^V^)) (2) 



where V(r) is the potential energy. For an isotropic glass under the assumption that 
all a.i are the same V e ^ becomes essentially a function of a scalar distance R, V e ^(R), 



and one finds 



1J a = § JdR g(R) V 2 (f3V eff (R)) and a free energy 



F(a) _ fideal 
Nk B T 



jideal _|_ 
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p f dR g(R) V e ff(R). Although a good approximation for large a, the self-consistent 
phonon theory is not accurate for low a. At low a, an expansion around the uniform fluid 
can be used, but this depends on knowing the direct correlation function which is not avail- 



able for network fluids. When no bonds are present, alternative 
first suggested by Tarazona|23( and generalized by Ashcroft 



y one can use the functional 
This functional takes the 



form N ^ T = f — J drp(r) [log(p(r)) — 1] + ^(r/) where $ is the interaction part of the 
free energy for a monatomic liquid or glass (which we take to be the Carnahan-Starling 
expression) and rj is the packing fraction (it pa 3 / 6). For the network forming liquids we 
will correct this free energy using a modification of the result of Erukhimovich|25|| . which 
provides for the entropy loss of forming bonds. 

For the network forming liquid and glass, we separate the Hamiltonian into bonding 
and non-bonding parts V = J2i,j bonded V b{ r ij) + £ fj - non-bonded Kifcj)- Each pair i and j 
is considered to be either bonded on non-bonded and thus have an interaction described 
by either V& or V^&. The molecular units which are considered to be bonded have a locally 
harmonic potential energy V&(ry) = Vo + ^-y b (r^ — d b ) 2 with d b the equilibrium bond distance. 
V n b, on the other hand, represents non-bonded interactions and therefore is taken to be a hard 
or soft sphere. Notice this potential does not explicitly describe, for any given configuration, 
whether a specific given pair should be considered bonded or not. This requires writing the 



26, 



. These manybody functions would 



potential using manybody switching functions 
reflect the saturation of chemical bonding tendencies (and would be needed for simulations.) 
For such a model the theory of associated liquids yields a (renormalized) equilibrium constant 
for the bonding reactions in terms of the underlying potential 28]. We will take the fraction 
of bonded pairs to be given and to be temperature and density independent. Of course, 
density and temperature dependent changes of the association equilibrium would be easy 



to include, but would distract from understanding how the glassy physics itself depends on 

temperature. Provided the association equilibrium constants are not paradoxically related 

to the density (i.e., decrease partially with p) the trends in any event will be as shown. 

For the network system, the self-consistent phonon calculation must distinguish the two 

types of interactions. For specificity we assume any given atom can form 0-4 bonds, in a 

tetrahedral pattern. Other spatial patterns are easily treated as well. At large a the effective 

potential coming from a bonded pair is essentially VJ, e ^(R) = -y + ^- (|R| — c^) 2 where the 

factor of 1/2 comes from the self-consistent formulation. We set the bond distance d\, equal 

to the hard sphere diameter for the non-bonded atoms and choose a zero of energy such that 

Vo = 0. In numerical calculations we used /?7 b = 300, approximately the value appropriate 
f~| 

for silica potentials [29]. The probability that a given atom has n-bonds, p n is given by 



Erukhimovich 25] p n = (^)r ra (l — r) 4_n with T = ^ and n b is the average number of bonds 
per atom. 

In our calculations the basic pair distribution for non-bonded interactions, g(R) is taken 
to be the Percus-Yevick radial distribution function modified for high density in the standard 
way. To account for the presence of bonded interactions that block the approach of other 
partners, g{R) for an atom with n bonds (denoted g n (R)) is scaled by max ^0, < ™^~ w j for 
distances less than the first minimum of g(R). < nn > p is the number of nearest neighbors 
at total density p, obtained by integrating the first peak of g(R). 

To each of the 5 types of atoms, those with 0, 1, 2, 3, or 4 bonds, we assign a mean a n . 
In this paper we use the uncoupled oscillator implementation of the self-consistent phonon 



theory 1J] to determine the free energy. Essentially, this approximation treats the glass as 
a collection of Einstein oscillators. The derivation in Stoessel and Wolynes's early paper is 
lengthy in detail but when applied to the present potential one obtains 



F scp /Nk B T = Y J p n V b ef n f (d b )+p I dRg n (R)V:i f (R) 

n=0 L J 

-In / rf r ^e- E "-«^ n -[ VV/3V/ ^ /(Rl)+Ia " b ]- ri (3) 

J V 

where V^J is the sum of V^^ for an atom with n bonds over its bonds (and includes 
correlations between bonds through the use of a tetrahedral lattice), V*[* is defined by 
Eqn. |21 using only the non-bonding potential energy and g n , and a™ b is the non-bonding 
contribution to a n . With a proper choice of local axes for each the argument of the 
exponential can be diagonalized, leading to 3 principal values of the vibrational tensor ai }Jl , 
a 2) n, and a 3jn . For the free energy to be consistent with Eqn. [H we must have a n satisfy 
3a' 1 = \ + ou^n + aj* Note that for a glass transition to occur, all 3 components of a n 
must be large. Given a bonding pattern, the {a} can be determined and, hence, the free 
energy. 

F scp is the free energy of the large a glassy minimum. We compute the free energy of 
the network liquid with a = using the Carnahan-Starling expression for a non-bonded 



liquid plus a modification of Erukhimovich's result for equilibrium network materials 
F liq /Nk B T = V(ti) + 4 [r In T + (1 - T) ln(l - T)} + In (f ) - f In (=*£) + Ik- whe re F h scp is 
obtained from Eqn. 01 with a = and only the bonding terms. It corresponds to the entropy 
loss from bonding computed by ErukhimovichjJ. 

The self-consistent phonon theory alone allows us to find pa, as the lowest density giv- 
ing a non-zero a, which thus corresponds to the onset of activated behavior while px, the 
"Kauzmann" density at which the free energies of the liquid and glass match and the con- 
figurational entropy vanishes also relies on the bonding correction to F(a = 0). The ratio 



Pa/ Pk dimensionlessly characterizes the thermodynamic aspects completely. To compare 
to the kinetic laboratory glass transition, we note the laboratory glass transition defined 
to occur when the viscosity reaches 10 14 p. Previous work using random first order theory 
predicts this to occur when the configurational entropy is about 1.0 ks- To make the trans- 
lation of our thermodynamic results to the laboratory transition density in what follows we 
will therefore mean by pc the density where the liquid and glass free energies differ by 1.0 
ksT per particle. The universality of configurational entropy at the laboratory transition is 
well confirmed, also by experiment. 

Pa, Pk, and pa as functions of the average number of bonds n b are shown in Fig. [TJ 
For rit > 3, p& vanishes, indicating that the rigidity percolation occurs when rib = 3 in 
our model. For comparison, Thorpe more precisely defining rigidity percolation as the 
absence of zero mode vibrational frequencies, finds rib = 2.4. For n& > 3, even at low 
densities the network is a rubber that can still vitrify at a higher densities. Fig. |21 shows the 
effective spring constants at these densities. For the hard sphere model|l4j the uncoupled 
oscillator approximation leads to a somewhat lower pa than simulation and to relatively 
small, but nevertheless finite, values of a a- Quantitatively, the more complete coupled 
oscillator approximation (which uses a chain summation to treat the glass as a collection 
of Debye oscillators |l^|) leads to a larger pa in better agreement with computer studies. 
So we would expect quantitative modification by taking better account of the coupling but 
expect the same trend, namely that as the amount of bonding is increased, the transition 
densities decrease and therefore the transition temperatures increase. This theory then 
explains why the addition of a non-bonding impurity into a networked glass will lower the 



glass transition temperature, as known since ancient times for silica 



In physical terms, 



the non-bonding transition densities decrease because, as seen in Eqn. as the contribution 
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from bonding to the vibrational force increases, the remaining contribution needed from 
non-bonding interactions for self-consistency decreases. 

It is clear that pa has a much sharper dependence on density than either p K or p G . 
The degree of bonding has a very direct effect on the vibrational force constant, being 
proportional to the average Laplacian of the effective potential, but only an indirect effect on 
entropy, entering through the more slowly varying (in space) Boltzmann factor. The theory 
explains the change of fragility with bonding. To see this we note that real liquids have 
both repulsive cores of finite strength and attractions which conspire in setting the equation 
of state. Yet assuming the attractions provide merely a smooth background potential the 
density dependences of this model can be translated into temperature dependences. For a 
soft potential V ~ r~ n we have the scaling relation T\jTi = (p2/pi) 3 ^ n - Using this relation 
for n = 12 yields Ta/T g and T K /T G as plotted in Fig . El Experimental measurements of 



31 



, while for the strong glass silica the 



Tk/Tq give 0.82 for the fragile glass orthoterphenyl 
ratio has recently been determined U to be 0.36. The uncoupled oscillator approximation 
results in a variation from 0.82 for the non-networked glass to 0.6 at the approximate rigidity 
percolation limit, a variation somewhat smaller in magnitude. Doubtless, the crudeness of 
this result is due to our use of a scaling approximation appropriate to a simple, purely 
repulsive model and better agreement with experiment would be obtained by including 
attractive forces in our estimates. We see that Ta/Tq increases dramatically with ri&, while 
Tr/Tq decreases modestly. Thus highly bonded liquids have large temperature ranges of 
nearly Arrhenius activated behavior while liquids with few crosslinks have smaller ranges of 
strongly non- Arrhenius activated behavior. 

For the uncoupled oscillator treatment a a falls rapidly with density to very small values 
as crosslinking increases. The universality of the Lindemann ratio previously assumed by 

9 



Wolynes and coworkers |3( would seem in doubt. Yet this is not so. First, we know the coupled 
oscillator approximation will give a larger pa therefore decreasing the range of variation of 
a a- For n b = 0, the coupled oscillator gives a value of a a ~ 100, while the uncoupled theory 



gives a a ~ 20|l4j. More important, even at the uncoupled level ock and «g are roughly 
constant with increasing bonding. A careful examination of the earlier analysis shows that 
these are the appropriate Lindemann-like ratios to determine the activation barriers in the 
strongly glassy regime. The predicted universal relation between D and AC P , relying on a 
universal «g, remains therefore sound. 

The combination of self-consistent phonon and liquid state approaches provides a micro- 
scopic picture of network glasses that explains the dependence of the characteristic glass 
densities and, less directly, temperatures on degree of crosslinking and this results agrees 
with experiment. Doubtless, many of the detailed liquid state approximations made here can 
be improved but it is unlikely such modifications will change the qualitative story. In partic- 
ular, the attractive forces that go beyond the soft-sphere model should probably be included 
to achieve quantitative agreement with experiment. This model combined with a theory 
of the association equilibria in forming the network provides a first principles treatment of 
glass properties starting from the underlying intermolecular forces. 

RWH acknowledges NSF Grant. No. 9977124 for support. 
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